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ABSTRACT 

We present the first identification of large-scale structures (LSS) at z < 1.1 
in the Cosmic Evolution Survey (COSMOS). The structures are identified from 
adaptive smoothing of galaxy counts in the pseudo-3d space (a, S,z) using the 
COSMOS photometric redshift catalog. The technique is tested on a simulation 
including galaxies distributed in model clusters and a field galaxy population - 
recovering structures on all scales from 1 to 20' without a priori assumptions for 
the structure size or density profile. Our procedure makes no a priori selection 
on galaxy spectral energy distribution (SED, for example the Red Sequence), 
enabling an unbiased investigation of environmental effects on galaxy evolution. 

The COSMOS photometric redshift catalog yields a sample of 1.5 x 10 5 galax- 
ies with redshift accuracy, Az WffM /(l + z) < 0.1 at z < 1.1 down to l AB < 25 
mag. Using this sample of galaxies, we identify 42 large-scale structures and clus- 
ters. Projected surface-density maps for the structures indicate multiple peaks 
and internal structure in many of the most massive LSS. The stellar masses 
(determined from the galactic SEDs) for the LSS range from M* ~ 10 11 up to 
~ 3 x 10 13 M . Five LSS have total stellar masses exceeding 10 13 M Q . (To- 
tal masses including non-stellar baryons and dark matter are expected to be 
~ 50 — > 100 times greater.) The derived mass function for the LSS is consis- 
tent (within the expected Poisson and cosmic variances) with those derived from 
optical and X-ray studies at lower redshift. 

To characterize structure evolution and for comparison with simulations, we 
compute a new statistic - the area filling factor as a function of the overdensity 
value compared to the mean at surface overdensity (/a(S/S(^)). The observa- 
tionally determined f^ has less than 1% of the surface area (in each redshift slice) 
with overdensities exceeding 10:1 and evolution to higher overdensities is seen at 
later epochs (lower z) - both characteristics are in good agreement with what 
we find using similar processing on the Millennium Simulation. Although similar 
variations in the filling factors as a function of overdensity and redshift are seen 
in the observations and simulations, we do find that the observed distributions 
reach higher overdensity than the simulation, perhaps indicating over-merging in 
the simulation. 

All of the LSS show a dramatic preference for earlier spectral energy distri- 
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bution (SED type) galaxies in the denser regions of the structures, independent 
of redshift. The SED types in the central 1 Mpc and 1-5 Mpc regions of each 
structure average about 1 SED type earlier than the mean type at the same red- 
shift, corresponding to a stellar population age difference of ~2 - 4 billion years 
at z = 0.3 to 1. 

We also investigate the evolution of key galactic properties - mass, luminosity, 
SED and star formation rate (SFR) - with redshift and environmental density 
as derived from overdensities in the full pseudo 3-d cube. Both the maturity 
of the stellar populations and the 'downsizing' of SF in galaxies vary strongly 
with redshift (epoch) and environment. For a very broad mass range (10 10 — 
— 1O 12 M ), we find that galaxies in dense environments tend to be older - this is 
not just restricted to the most massive galaxies. And in low density environments, 
the most massive galaxies appear to have also been formed very early (z > 2), 
compared to the lower mass galaxies there. Over the range z < 1.1, we do not 
see evolution in the mass of galaxies by more than a factor ~ 2 separating active 
and inactive star-forming galaxy populations. 

Subject headings: cosmology: observations — cosmology: large-scale structure of 
universe — cosmology: dark matter — surveys 



Introduction 



The Cosmic Evolution Survey (COSMOS) is intended to probe the evolution of galaxies, 
AGN and dark matter in the context of their cosmic environment (large-scale structure - 
LSS). The survey area samples scales of LSS out t o ~ 50 ~ 100 Mpc at z > 0-5. This 
corresponds to 9 times the area of GEMS and EGS (IRix et al.ll2004i iDavis et al.ll2006l ). the 
next largest HST imaging surveys. The COSMOS area (jScovilld l2006al ) is expected to have 
a 50% probability of having one 10 14 M halo (dark and luminous matte r) within every Az 
~0.1atz~l-2 (based on ACDM simulations; see lBenson et al.l (120011 )): lower mass halos 
(~ 10 13 M ) are ~ 20 times more abundant and therefore will be seen in every Az ~ 0.1. 
A major theme for COSMOS is the effect of cosmological environment on the evolution of 
galaxies and AGN. The identification and measurement of LSS are therefore a prerequisite 
to many aspects of science with COSMOS since the large-scale structures define the local 
environment. In COSMOS, the local galaxy number and mass densities can be compared 
with the total mass densities determined from weak lensing tomography and ho t X-ray 



emitting gas in the virialized parts of LSS having clusters of galaxies (see below, iMassey 
20061 : iFinoguenov et alJliooS ). 
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The identification of LSS from the observed surface-density of galaxies requires some 
means of discriminating galaxies at different distances along the line of sight; otherwise, the 
increased shot noise in the galaxy counts reduces the signal-to-noise ratio for the large-scale 
structure. The better the redshift or distance discrimination, the easier it is to see low density, 
large-scale structures. (This is true down to the point that the internal velocity dispersion 
of the structure is resolved.) For LSS finding, line-of-sight discrimination is usually accom- 
plished using: 1) color selection (e.g. using broadband colors to select red sequence galaxies, 
Gladders fc Yed 120051 ); 2) photometric redshifts based on the bro adband spectral energy 



distribution (SEP), or 3) spectroscopic redsh i fts (see Appendix |A~1 Ivan de Weygaertl Il994 



Postman et al. 1996 ; Schuecker fc Boehringer [199J ; iMarinoni et all 120021 ) . Color selection 
is not used here since it would preclude investigati on of correlations b e tween environmen - 
tal density and galaxy SED or morphological type (IDressler et al .1119971 : ISmith et al.l 120051 ). 
With color selection, the defined large-scale structures would be a priori biased to a particular 
galactic SED type (e.g. early type galaxies if the red-sequence method is used). High density 
structures may well be rich in red galaxies of early morphological type, but exploring the de- 
pendence of galaxy type on environmental density requires that the environment be defined 
or identified without bias toward specific galaxy types. Ultimately, the use of spectroscopic 
redshifts to determine distances is more desirable, provided a sufficiently large sample exists 



(iLe Fevre et al.ll2005t iGerke et al.ll2005t Meneux et al.ll2006l : ICooper et al.ll2006l : ICoil et al. 



20061 ). Th e z-COSMOS spectrosco pic survey will yield ~ 30,000 galaxies when completed 



in ~2008 (iLillyl 120061 ; llmpeyi 120061 ); however, at this time the spectroscopy is much more 
limited and we must rely on the alternative approaches for line-of-sight discrimination. 

In this paper, we identif y LSS in the 2D°COS MOS field using the extensive COSMOS 
photometric redshift catalog (IMobasher et al.ll2006l ) to analyze the galaxy surface density in 
redshift slices out to z = 1.2. The galaxy samples used for this work and their completeness 
are discussed in Section [2J We use an adaptive smoothing technique (see Appendix |A| to 
identify areas of significantly enhanced galaxy surface-density. For each significant peak 
in the smoothed surface-density pseudo 3-d cube, we find all connected pixels to delineate 
a sample of 42 structures (Section [3]). For each structure, we estimate the dimensions, 
number of galaxies, and mass (from the broadband fluxes of the galaxies). The stellar mass 
distribution of identified COSMOS LSS extends from 10 11 up to 3 x 10 13 M . The relative 
amount of structure at different overdensities is analysed in Section H] and compared with 
results from the Millennium Simulation. 

We then investigate the evolution of galaxies with respect to their location in the 42 LSS 
in Section and environmental density in Section EJ Strong variation of the SED type and 
star formation activity is seen with both redshift and environmental density. The maturity 
of the stellar populations and the 'downsizing' of SF in galaxies is also strongly varying with 
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epoch and environment (Section l6.4.3l) . (Adopted cosmological parameters, used throughout, 
are: H = 70 km s' 1 Mpc" 1 , Q M = 0.3 and Q A = 0.7) 



2. Photometric Redshifts and Sample Selection 



For this investigation we use photometric redshifts to separate the galaxy population 
along the line-of-sight. It is vital for the analysis that the binning in redshift be matched 
to the accuracy of the redshifts. Using binning that is finer than the redshift uncertainties 
distributes the galaxies in a single structure over multiple redshift slices and thus reduces 
the signal in each slice. Conversely, bins of width larger than the redshift uncertainties will 
increase the shot noise associated with the background surface-density of galaxies, relative 
to the large-scale structure signal. 

The COSMOS photometry catalogs were generated from deep ground-based optical 



imaging at Subaru (jTaniguchi et al.l 120061 ) and CFHT; they are com bined with shallower 
near infrared imaging from NOAO (KPNO and CTIO) fjCapakl 120061). The resulting pho - 



jasher et al.l 120061 ) . 
Lilly! tod )), there 



tometric redshift catalog contains 1.2 million objects at Iab <26 (IMo 
For approximately 900 objects (all with Iab < 24 mag from zCOSMOS 
exist spectroscopic redshifts; after removing 'catastrophic' outliers (~1% of the sample), 
the offsets between the photometric and spect roscopic redshifts have o z j{\ + z) ~ 0.03 or 
Azfwhm / 0- + z) < 0.1 (IMobasher et al.l 120061 ) . Since the spectroscopic - photometric red- 
shift comparison is limited to a small sample of mostly brighter objects, we will instead 
use the 'goodness' of fit from the photometric redshift determination for a more general 
assessment of the redshift accuracies. This represents an internal dispersion and hence is 
likely to be lower than the true uncertainty which includes systematic effects; nevertheless, 
it does provide a characterization of the dependence on redshift and magnitude cutoff for 
any sample selection. 

T he photometric redshifts were derived using the Bayesian Photometric Redshift method 
(BPZ: iBenson et al.ll2000l ; IMobasher et al.ll2006l ). The fitting presumes 6 basic spectral en- 
ergy distributions (SEDs) and for the photometric redshift catalog used here, there was no 
assumed 'prior distribution' for the galaxy magnitudes. The dust extinction within each 
galaxy was also a free parameter in the redshift fitting. For SED types to 2, a Galac- 
tic extinction law was assumed: for SED types greater than 2, a Calzetti law was used 
(IMobasher et al.l 120061 ) . The fitting outputs the most probable redshift, 68 and 95% con- 
fidence intervals, the intrinsic SED type and the absolute magnitude (My). In Figure [H 
the redshift distribution and uncertainty in the redshift fits are shown as a function of red- 
shift and i-band magnitude cutoff in the sample. The SED classifications used here are : 1 
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(E/SO), 2 (Sa/Sb), 3(Sc), 4(Im), and 5-6(starburst). The quantity A z(50%)/(l+z) shown 
in Figure [lb is the mean value (at each z) of the width in z containing 50% of the probability 
distribution. (This full- width is 1.3 x a for the derived fits, assuming a Gaussian uncertainty 
distribution.) Also shown is the photometric redshift uncertainty as a function of apparent 
magnitude cutoff - based on these curves, we adopt a cutoff Iab = 25 mag. The uncertainties 
plotted in Figure [1] indicate that the bin width for identifying large-scale structures should 
be approximately Az ~ 0.1 — 0.15 up to z ~ 1.2 and Az ~ 0.25 up to z ~ 2.5 for the chosen 
magnitude limit. 

Throughout the investigation here, we use the galaxy rest frame SEDs to characterize 



the galaxy type, rather than the observed morphologies. ICapak et al.l (120061 ) find a tight 
correlation between SED and morphology as measured by the Gini and Compactness mea- 
sures. We have also correlated the Sersic indices measured using GALFIT for 5,000 bright 
(I < 22 mag) galaxies in COSMOS at z = 0.2 to 0.8. The SED type 1 (E/S0) is strongly 
correlated with an R 1 / 4 law (Sersic n = 4); however, there exists a large dispersion in the 
Sersic indices for the later SED types. This large dispersion probably reflects both the real 
dispersion in bulge to disk ratios of later galaxy types and the difficulties of measuring the 
morphologies for faint galaxies at high redshift. Throughout this investigation we will use 
the SED types, derived from the photometric redshift fitting, to classify the galaxies since 
the SEDs are more readily classified for faint galaxies than the morphology. 

The COSMOS photometric redshift catalog also includes galactic stellar masses, derive d 



from the absolute magnitude, SED type and redshift of each galaxy flMobasher et al.ll2006l ). 
Approximate estimates for the star formation rates (SFR) were also derived using the SED 
type, absolute magnitude and fitted extinction. The intri nsic UV continu um (corrected for 



extinction at A = 1500 A) was used to estimate the SFR (jMobasherl 120061 ) . 



2.1. Galaxy Samples 

Although the COSMOS photometric redshift catalog contains over a million objects, we 
impose selection criteria to yield a more reliable galaxy sample for structure identification 
- i.e. galaxies with the best photometric redshifts, detected in several bands, and with 
significant intrinsic luminosity. We thus restrict the analysis to : 

25 mag > Iab > 19 mag (la) 

My < -18 mag (lb) 

We also require that each object be detected in at least 4 bands and the SExtractor stellarity 
parameter be less than 0.95. The former (in addition to the 25 mag cutoff - Eq (la)) limits 
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the sample to galaxies with accurate photometry; the latter removes objects which are likely 
to be stars or QSOs. Similarly very bright objects are also excluded by condition (la) since 
they are likely to be stars. Since the fraction of galaxies with dominant AGN is probably 
not large, their exclusion should not significantly affect the large-scale structure definition. 
Condition (lb) removes galaxies which have low absolute luminosities and presumably low 
mass. These criteria yield the galaxy samples listed in Table [1] along with the adopted 
redshift binning for the LSS identifications presented below. In Table HJ we also provide a 
breakdown of the samples with respect to SED type from the photometric-redshift fitting. 
The majority of the analysis in this paper refers to the low-z sample in Table [H 



2.2. Galaxy Selection and Completeness with Redshift 



At larger redshifts, the galaxy sample used to define LSS will be incomplete at low 
luminosities (masses) and we evaluate here the severity of this effect with two approaches : 
evaluating the cutoff L* (the characteristic luminosity at the knee in the Schecter luminosity 
function) as a function of redshift and comparing the galaxy mass functions as a function of 
redshift for the sample described in Section 12.11 

Figure [2] shows the distributions of absolute magnitude (left) and stellar masses (right) 
for galaxies in our sample as a function of redshift. On the left panel, the lines indicate 
the expected absolute magnitudes for L* and L*/10 galaxies assuming a typical, passive 
evolution brightening of 1.2 mag from z = to 1.2. Here one sees that the sample easily 
goes down to L*/10 out to z = 1.2 with Iab < 25 mag (the lower envelope in the grayscale). 
Incompleteness does set in at less than L*/10, but typically less than 30% of the total 
luminosity is contained in these galaxies for a Schecter luminosity function. 

Alternatively, one can compare the distribution functions of galactic stellar masses as 
a function of redshift to assess the incompleteness (assuming, to first o rder, that th e mass 
function is not strongly varying over this redshift interval - e.g. < 2x; see lBorch et al.l (120061 ) 
and below). The derived mass functions for galaxies entering the sample (Section |2.ip used 
here for LSS definition are shown on the right panel of Figure [2j These mass functions are 
in agreement in both s hape and absolute value with previously determined m ass functions 
for this redshift ranee (IDrorv et al.ll2005l : iBorch et al.ll2006l ; iBundv et al.ll2005T ). The higher 
noise seen in the low z mass functions is due to the much smaller volume and hence smaller 
number of galaxies sampled. 



The total number of galaxies and total mass of galaxies per unit comoving volume 
were evaluated by integrating the distribution functions shown in right panel of Figure [2] 
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at masses above 10 9 M & . In columns 2 & 3 of Table [21 the totals are divided by the 
comoving volume in order to assess the count and mass incompleteness relative to this 
maximum redshift bin (see Table [2]). The falloff seen in the mass functions for z > 0.7 at 
^-stellar < 5 x 1O 9 M is probably due to incompleteness at the apparent magnitude limit 
Iab < 25 mag for our sam ple. Incompleteness at this magnitude limit is quantified for 
COSMOS in IScovilld (j2006al ). In terms of integrated stellar mass for galaxies above 1O 9 M , 
our sample is at most missing only 10% (Table [2]) of the total mass, relative to the lower 
redshift bins which are more complete (e.g. z ~ 0.5). In the analysis below we will not 
correct for this incompleteness unless noted explicitly, given the fact that it is probably not 
large and the uncertain assumptions of constancy in either the mass- or luminosity- function 
which would be required. 



2.3. Pseudo 3-D Surface-Density and Noise Estimates 

The adaptive smoothing algorithm we employ here is designed to analyze redshift slices, 
each of which represents the surface-density of galaxies (£) in a redshift bin. The custom- 
built algorithm is formally described along with test results in Appendix |A] As noted in 
Section[2l the width of these slices will be Az = 0.1 and 0.25 for the low and high-z samples 
respectively. However, given that the different galaxies may have quite different widths for 
the fitted redshift probability distribution, insertion of each galaxy into the 3-d cube (a,5,z) 
as a delta function (S(z — z p k)) at the most likely redshift would not optimally weight the 
galaxies with the most accurate photometric redshifts. Instead, we populate the 3-d cube 
with a Gaussian distribution in z for each galaxy. The Gaussian dispersion was taken from 
the high and low-z 68%-confidence limits from the photometric redshift fit - specifically 
(?z = (ztiigh 68% — z iow 68%)/2). Thus, in the adaptive smoothing procedure, galaxies which 
have a large uncertainty in their derived redshifts will have relatively low weight, because 
they will be spread over a larger range in the redshift dimension. And galaxies with tighter 
redshift fits will be treated more significantly. One concern might be that this tends to 
prefer structures defined by early type galaxies which have a strong Balmer break and thus 
small redshift uncertainty. As a test, we also used the adaptive smoothing on a 3-d cube 
with the galaxies located as points (rather than a probability distribution) at their most 
probable redshift, Since this test yielded structures similar to those shown here, we prefer 
employ the probability distributions to take account of the redshift uncertainties. The cube 
being analysed is therefore the 3-d 'probability' surface- density of galaxies, not the galaxies 
as discrete points in 3-d space. For the adaptive smoothing, the required noise estimate (a) 
is taken as the counting uncertainty, i.e. the square root of the galaxy surface-density cube. 
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The square COSMOS field is 1.4xl.4° in s ize; the comoving volume out to z = 1.1 in 
the low-z sample is ~ 10 7 Mpc 3 (jScovilld l2006bl ) . For the adaptive smoothing, we use a grid 
of 300x300 in (a, 5). The angular resolution in the smoothing is therefore ~ 17". A typical 
redshift slice with Az =0.1 contains ~ 1-3 x 10 4 galaxies (see Figured]); each cell will 
therefore be populated by ~ 0.2 galaxies, on average. Significantly higher computational 
resolution is therefore not warranted. 



3. Adaptively Smoothed Surface-Density 

The galaxy surface- density derived using the procedure described above and in Appendix 
lAlis shown in Figure [3] for redshift slices with Az = 0.1, spaced by Az = 0.1 for the Low-z 
sample of galaxies. In the last two panels of Figure [3] we show two of the higher redshift 
slices with Az = 0.25. We leave further analysis of the high redshift g alaxies to a later paper 



since deeper near infrared and Spitzer IRAC imaging (jSandersI 120061 ) is required for higher 
accuracy photometric redshifts. 

The surface- density plots show a large number of very significant large features - espe- 
cially at z = 0.35, 0.75 and 0.85. And at every redshift numerous small groupings of galaxies 
are seen. There is a definite trend towards increasing complexity of structure (dumpiness) 
at higher z. This is to be expected since structures at high z (earlier epochs) are dynamically 
younger and expected to be less relaxed. 



3.1. Structure Identification 



From the derived surface density in the 3-D cube (a, 5 and z), we define preliminary 
LSS starting from > 1 0a peaks, finding all c onnected pixels above the la noise. Using an al- 
gorithm developed by lWilliams et all (119941 ). approximately 140 local maxima are identified 
and their connected pixels catalogued. When multiple local maxima are found in proximity, 
the neighboring pixels are associated with the nearest local peak - this can result in subdivi- 
sion of structures which have multiple peaks (real or noise). The maps of the 140 preliminary 
structures are therefore checked for possible recombination into composite structures. The 
decision to recombine was based on : whether the individual components were touching in 3- 
d space; their borders meshed; and their proximity in the 3-d space was unlikely by chance. 
This is somewhat subjective but a more physically justified recombination would require 
spectroscopic redshifts with accuracy similar to the virial v elocities of the groups. This will 
become possible with the COSMOS spectroscopic surveys (jLillyl 120061 : llmpeyl 120061 ). 
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Forty-two recombined, independent LSS were found from the procedure described above 
(i.e. identification of all local maxima and the subsequent recombination). The surface den- 
sity of each structure, integrated in the z-dimension, is shown in Figure |U Many of the most 
populated structures show complex structure with multiple peaks in a and 5. The structures 
are ordered in terms of decreasing number of 3-D pixels so the most complex and massive 
structures are those with lowest LSS number. In Table [3] measurements for the structures are 
tabulated, including the location of peak density a and S, centroid redshift, sizes, number of 
galaxies and mass. Figure shows the projection of all LSS on to the plane of the sky, i.e. 
integrating in redshift; a finding chart for the LSS within the COSMOS field is provided in 
Figure [61 The galaxy stellar masses are obtained from the photometric redshift fit, which 
yields an absolute magnitu de and an SEP type f or the galaxy from which a stellar mass-to- 
light ratio can be inferred (IMobasher et al.ll2006l ). An important consideration in measuring 
the structure parameters is possible contamination from the background galaxy population. 
This contamination is estimated from the mean surface- density in each redshift slice. Then 
for each galaxy, seen in projection within the area of a given LSS, the probability that it is in 
fact associated with the structure is given by the ratio ^LSs/C^LSS + ^background)- In calculat- 
ing the cluster mass, the mass of each galaxy within it is multiplied by this local probability. 
Thus, the derived masses are corrected for foreground/background contamination. In the 
first column of Table [31 we list ( in par enthesis) the ID numbers of possibly associated X-ray 
clusters from lFinoguenov et al.l (120061 ). The surface densities at the cores of these structures 
are simil ar to those seen in the Rich Cluster and Subaru/XMM Deep su rveys (S ~ 10 20 



Mpc" 2 ; lO'Helv et al.lll998l : iPimbblet et al.ll2006l : iKodama et al.l 120041 ) but their exten ts at 



lower density go out to > 10 Mpc. Spectroscopic redshifts from zCOSMOS (ILillyi 120061 ) will 
be needed to confirm the coherence of these more extended structures. 



3.2. Radial Profiles : Clusters vs. Structures 

Without spectroscopic redshifts it is impossible to determine which structures are in fact 
gravitationally bound. However, the spatial distributions of galaxies within the structures 
suggest that many of the LSS are 'relaxed' clusters. In Figure [3, we plot the azimuthally 
averaged projected radial distribution of galaxies in each structure. For each LSS, radii 
were calculated from the position of peak number density. In several of the structures with 
significant secondary peaks the radial structure is not monotonically decreasing (e.g. # 10, 
13 and 18). Since the largest, most complex structures have the lowest LSS numbers, these 
are most likely large-scale structures with multiple clusters (LSS # 1-8). All of these are 
at the higher redshifts; this is due to the greater volume sampled at z > 0.5. They may 
eventually relax to form a centrally-concentrated cluster. Conversely, LSS # 30 - 42 all 
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appear fairly symmetric in their radial distributions and with size 1-2 Mpc, similar to 
those of present day galaxy groups (< 20 members) or small clusters. 

Most of the structures can be fit by a power law surface density of roughly r _1 , within 
the central few Mpc, implying that the physical density is ~ r~ 2 - similar to what is usually 
found for local clusters such as Perseus. There are however some LSS where the density 
dependence steepens in the central regions; this could reflect the presence of an unrelaxed, 
outer 'infall' region. A better understanding of the nature of these density profiles and their 
variations will require better kinematics from spectroscopic redshifts. 



3.3. Richness 



The last column of Table [3] provides an esti mate of the R ichness of the structures using a 
measure similar to that used for galaxy clusters (lAbelllll958l ). For each structure, the central 
surface density of galaxies at R < 1 Mpc is listed with the background number counts 
subtracted off. Radius is measured from the location of peak surface density as in Figure [7J 
We first find the 3rd brightest (M3y) cluster galaxy and then count all galaxies brighter than 
M3y + 2 mag. (M3y is always brighter than -20 so this procedure is not conflicting with 
the cutoff in Equat i on lb .) The standard procedure for local clusters employs R < 1.5 Mpc, 
but iPostman et al.l (119961 ) find that for a typical cluster profile N(i? < 1 Mpc) / N(i? < 1.5 
Mpc) ~ 0.72 ± 0.05 , so the estimates given in Table [3] can be scaled up by approximately 
1.39 to make them comparable to the standard Abell Richness criteria. The distribution of 
Richness parameters is shown in the right panel of Figure [HI Most of the cores of COSMOS 
LSS fall in Richness classes 1-3. 



3.4. Structure and Cluster Masses 



Figure [8]shows the distribution of photometrically derived stellar masses (IMobasher et al 

20061 ) for the LSS between 5 x 10 11 and ~ 3 x 10 13 M . This distribution is clearly subject 



to significant Poisson and cosmic variances as discussed below (Section l3.5p . The observed 
distribution increases toward low mass, but there are 4 with masses exceeding 10 13 M & . On 
the high mass end of the da rk matter halo mass spectrum, the expected number distribution 
is N{m)dm oc m~ ls (e.g. iBenson et al.l (120011 ) ). The distribution of total stellar masses 



shown in Figure [8] is much less steep, but at this point the ratio of stellar to dark matter 
mass as a function of halo mass and z is not known. 



The highest mass structure is LSS #1 with a stellar mass of 2.3x1 13 M ; clearly 
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this is a super-massive structure, equivalent to that of the Coma cluster if allowance is 
made for a dark matter contribution. The mass in LSS #1 at z ~ 0.74 is distributed over 
scales ~ 10 Mpc. In fac t , the structure appears to be aggregating around a central cluster 
( IGuzzo Cassatal 120061 ; ICassatal 120061 ) and is therefore possibly forming a super-massive 
cluster like Coma. LSS #1 is also detected i n the weak lensing shear analysis (jMasseyl 120061 ) 
and in the X-rays (jFinoguenov et al.l 120061). LSS # 17 seems to exhibit a very complex 
sub-structure as discussed in detail by (jSmolciq 120061 ) . Within the inner ~ 2 Mpc of LSS 
#17 there are at least 4 X-ray luminous clusters and one X-ray quiet overdensit y at the 
same redshift . One of the clusters hosts a wide-angle tail (WAT) radio galaxy which ISmolcic 
( 120061 ) discuss as a tracer for assembly of this complex cluster. They argue that the structure 
is in the process of formation and estimate that the mass of the final cluster, after merging 
of all sub-components, will be ~ 20% of the Coma cluster mass. 

The LSS masses listed in Table [3] are for just the stellar masses as derived from the 
observed galaxy fluxes using a m ass-to- light ratio , base d on the best fit SED from the pho- 
tometric redshift determination (IMobasher et al.ll2006l ). The total masses, including non- 
stellar or non-luminous baryons and dark matter, are at least an orde r of magnitude greate r 
for n R = 0.025/r 2 and Q M = 0.3 with H = 70, Vt M /Vt B = 6.1 fcolb fc Turr^rl ll99oh . 



Hoekstra et al.l (120061 ) analysed the weak-lensing maps for a sample of individual galaxies at 
0.2 < z < 0.4 to estimate the total virial masses and baryon fraction in the stars. They find 
a virial-to-stellar mass ratio M V i r /M* = 20 — > 40, depending on the assu med stellar IMF 
and M*/M hnrynn = 14% (early type galaxies) to 33% (late type ga laxies) (Hoekst ra et al. 



20061 ). Similar results were found for lower redshift SPSS galaxies bvlGuzik k, Seliakl ( 2002) 



and i n semi-analytic simulations ( Kauffmann et al. 1999 ; Bemtez 200ol ; van den Bosch et al. 



20031 ). Since th e mass-to-light ratio (including dark matter) is found to increase as ~L 15 
at low redshift (IHoekstra et al.l 120061 ) and presumably, the stellar baryon fraction is lower 
at higher redshift, we adopt M vir /M* = 50 as a reasonable lower limit for the LSS listed 
in Table [3] and a more likely yalue might b e ~ 100. For LSS #1 which has been ana lysed 
in detai 
opti cal ( 



using weak lensing (|Massevl 120061 ) , X-ray emission (jFinoguenov et al.l 120061 ) and 
Guzzo Cassatal 120061) . the apparent ratio of total mass to stellar mass is ~ 50 to 



100 (IGuzzo fc Cassatal 120061 ). The total masses for the LSS are therefore likely to be in the 
range 10 13 to ~ 3 x 10 15 M P 



I©. 



3.5. Variances in Distributions 



The mass distributions derived for the LSS are subject to both shot noise, due the small 
number of structures within each redshift slice, and to the cosmic variance that characterizes 
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the mass distribution on very large scales. To estimate the resulting unc e rtaint ies in our 
LSS mass distributions we follow the method described by ISomerville et al.l (120041 ) . We first 
calculate the volume and total mass (M(vol)) contained in each redshift slice (an area of 2.5 
deg 2 for the photometric redshift cata log used here). The co smic variance on this scale is 
given by am {vol) (see Figure 3-right in ISomerville et al.l (120041 )) and the relative variance in 
number counts for halos of mass in a redshift slice is then : 



a{M h ,z) = {a\ 



shot) 



,1/2 



a(M h ,z) = {{b{M h ,z) a M (voi)) 2 + l/N h f' 2 



(2a) 
(2b) 



where b(Mh, z) is the bias for the halos of mass Mh at redshift z, calculated as inlSheth fc Tormen 



mm, 



1 9 9sj) ■ and is the average number of halos in the slice. As explained in Somerville et al. 
, (Jm(voI) should be a slight overestimate of the true cosmic variance, since the volume 



in each slice is much deeper than it is wide, and thus along its z-axis the slice samples much 
larger scales where the variance is smaller. 



As suggested by lMo fc White! (119961 ): ISomerville et al.l (12004J ) . we can identify the appro- 
priate mass range for halos corresponding to the LSS by integrating the halo mass function, 
normalized to the survey volume in the redshift slice, down to a threshold mass which yields 
a total number of halos matching the observed number of LSS. This assumes that the de- 
tected LSS corresponds to the most massive halos on a roughly one-to-one basis, and thus 
that the overall abundance of LSS indicates the characteristic mass-scale of halos with which 
they are associated. 

In Table HJ we summarize the expected relative variances (shot, cosmic and total) in 
each slice for two different mass ranges of DM structures (10 13 - 10 14 and 10 14 - 10 15 M ). 
These variances were calculated for the co smological paramete rs specified in Section [TJ and 
with as = 0.74 (the WMAP 3-year value, ISpergel et al.l (120061 )). Comparing the expected 
numbers of halos for the two mass ranges with the observed numbers of LSS, it is most 
reasonable to identify the observed LSS with the higher mass range halos, i.e. 10 14 - 10 15 
Mq, for which the expected number is ~ 37. (This identification is only approximate since 
clearly some of the observed structures are much less massive.) 

Based on the results shown in Table HI we expect that the shot or Poisson noise and 
the cosmic variance are quite comparable for the mass range of the LSS sampled here at all 
redshifts. The total combined relative variance is expected to be in the range 0.4 to 0.6; 
i.e., for the very small number of very high mass structures the derived mass and number 
distributions will have typical uncertainties of ~50% for each redshift bin. 
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3.6. Redshift distributions of LSS 

Redshift distributions of the structures in two ranges of LSS stellar mass, M* are shown 
in Figure [9] As discussed in Section 13.51 uncertainties due to Poisson and cosmic variance 
are comparable and the expected total relative variance in these number distributions is 
40-60% (i.e. <Jn/N ~ 0.5). Also shown is the relative comoving volume (dotted line) for the 
redshift bins. The redshift distribution of total mass (EMlss) within structures with stellar 
masses in the range M* = 10 12 — > 10 13 M© is similar (within the expected 50% variance) to 
the dotted curve showing the variation of comoving volume sampled. This suggests that we 
are recovering structures in this mass-range without a strong redshift-dependent selection 
bias. The higher mass LSS (M, > 10 13 M ) exhibit an apparently steeper falloff at low z 
but this is not statistically significant given the small volume sampled. These results are 
consistent with a lack of dramatic evolution in the overall mass fraction for the most masssive 
struc tures out to z = 1 (a lookback time of ~ 8 Gyr) - as is also seen in ACDM simulations 
(e.g. iBenson et al.l (120011 )). However, this result is certainly not strongly constraining given 



the large variances. 

The derived mass function for the COSMOS LSS can be compared with previously 
derived mass distributions for galaxy clusters, mostly for clusters at lower redshift. The 
cumulative mass function (n(> M) = f™n(m) dm) is shown in Figure [TU] for the 42 COS- 
MOS LSS. The error bars are taken from the Poisson noise in each bin of width 2 x 1O 12 M . 
The expected cosmic variance is comparable to the Poisson noise (see Section 13.51 and Ta- 
ble SJ); it is not explicitly included here since it is dependent on the adopted cosmo logical 
parameters (in particular er 8 ) and on correct identification of the associated DM halo mass 
range. The full error bars are likely ~ v(2) — 40% bigger than shown when including 
cosmic variance. For the COSMOS sample volume we adopt 1.5 xlO 7 Mpc 3 out to z = 1.1. 
Als o shown are the mass functi ons derived from optical and X-ray studies, as summarized 



by (IReiprich fc Bohringer! 120021 ) . The mass function for masses within the Abell radius of 



each cluster is M^; masses within the regions with density exceedin g 200/? r and 500p r are 



M2nr) a nd M.^o (where p r . is t he cri tical density for the universe, see IReiprich fc Bohringer 



20021 ). IReiprich Sz Bohringer! (120021 ) derive the total mass including the dark matter and we 
have scaled their masses down by a factor of 100, i.e. assuming a stellar mass fraction of 
1% of the total mass (baryons plus dark matter). We have also scaled to h = 70 (used here 
throughout) from their h = 50. 

Figure [10] shows reasonably consistent number densities (pe r co-moving volume) betwee n 



the COSMOS LSS and the previous studies as summarized by (IReiprich fc Bohringer! 120021 ). 
given the somewhat uncertain ratio of total to stellar masses (taken as 100 for Figure [TO]) . 
It should be noted that the mass function within the Abell radii (labelled M^ in Figure [TU]) 
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also closely approximates the mas s function derived by (jBahcall &: Cenlll993l ). As noted in 
Section 13. 4[ [Hoekstra et al.1 (120061 ) determined a value of up to 40 for this mass ratio based 
on lensing measure for clusters at z = 0.2 to 0.4. The somewhat higher value, found here in 
order to achieve agreement in the local mass function shape, might indicate that the fraction 
of baryons in stars is less at the higher redshifts (z ~ 0.2 to 1.1) sampled in the COSMOS 
LSS. Alternatively, the COSMOS LSS measurem ents refer to more ext ended, lower density 
structures and filaments than those sampled by iHoekstra et al.1 (120061 ) and the conversion 
efficiency of baryons into star is very likely dependent on environment in the context of 
ACDM models. 



4. Comparison of Structures with ACDM Simulations 

ACDM simulations provide quite specific and relatively confident predictions for the 
growth structure in the dark matter (DM) as a function of redshift, given a specified set 
of cosmological parameters. On the other hand, the formation and evolution of the visible 
galaxies within the DM structures has relied on semi-analytic models or prescriptions for star 
and AGN formation, stellar evolution and feedback processes. These semi-analytic models 
and the predicted distributions of galaxies in ACDM have been mostly constrained from low 
redshift galaxy surveys. Relatively little constraint or testing of the semi-analytics vis-a-vis 
the DM LSS has been done at high redshift (i.e. z > 0.2). In this section, we compare in 
some detail the distributions of galaxy overdensities seen in the COSMOS field out to z = 
1.1 with those predicted in simulations. 

In particular, we will compare the relative volumes (or areas) occupied by observed 
structures of overdensity with the simulation predictions - as a function of redshift. As 
time progresses, the fraction of volume with high overdensity will increase and the max- 
imum overdensity should increase at lower redshift. This measure of structure evolution 
enables significant comparison between the simulations and the observed universe, avoiding 
the azimuthal averaging which is inherent is a correlation function analysis. The structures 
are expected to filamentary and therefore are not circularly symmetric; they may also have 
multiple characteristic scales. For the same reasons, we have employed the adaptive smooth 



ing technique developed here rather than matched filter algorithms (IPostman et al.1 11996 



Schuecker fc Boehringerill998l ) which are obviously well-adapted to the central, high density 



core structures but less appropriate to extended filamentary s t ructu res. Angular correlation 



functions for the COSMOS field are presented in iMcCrackenl (120061 ) 



Figure [TT] shows the fractional cumulative area with galaxy surface density greater than 
£/ < X > where < X > is the average in each redshift slice. The three colored curves show 
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overdensity filling factors for the redshift ranges 0.2 to 0.5, 0.5 to 1.1 and 0.2 to 1.1. These 
curves were computed from the adaptively smoothed overdensities shown in Figure [3] divided 
by the mean background surface density (E B) given in the top legend for each redshift slice). 
The level of the background is dependent on the density of true 'field' galaxies and on the 
redshift accuracy - thus in comparing with simulation predictions below, we also convolve 
the redshifts of galaxies in the simulation with a Gaussian of z-width matched to that of 
the observational photometric redshifts. Figure [TT1 exhibits the basic characteristic expected 
for structure growth as a function of redshift - higher overdensities occuring at later times 
(lower z) and a larger fraction of the area in overdense regions as time progresses. 

A quantitative comparison can be made with the Millennium Simulation. Mock catalogs 
were constructed using the Virgo Consortium's Millennium Simulation and the Galform semi- 
analytic model of galaxy formation. Dark mat ter and merger tre es were extracted from the 



Millennium Simulation using the techniques of iHelly et al.l (120031 ) , utilizing all halos of > 20 



particles. These merger trees are f ed through the Galform semi-analytic model (using the 



parameter set of iBower et al.l (120061 ) to populate the simulation with galaxies at all redshifts. 
We did not have access to proprietary lightcone data from the Millennium Simulation, so 
the mock catalogs were constructed taking cubes from the Millennium Simulation at z=0.3, 
0.5, 0.7, 0.9 and 1.1. Regions with two sides equivalent to 1.4° and one side extending 500 
Mpc/h were extracted at each redshift. Galaxies were selected to have My < —18 mag and 
19 < i <25. The mock based on the Millennium- Virgo semi-analytic (black curve in Figure 
[TTT) is in remarkably good agreement with the mean curve determined from the observations 
for z = 0.2 to 1.1 (solid green line). Jackknife tests were done, splitting the data in half 
and the variances are typically < 20% for most values of the overdensity - this provides 
a limited estimate of the uncertainties. The overall area filling factors in the observations 
and theory track each other within a factor of ~ 2. It does appear that the theoretical 
curve does not reach as high overdensities as the observational curve - possibly indicating a 
significant discrepancy on small scales. Since what is being measured in both the observed 
and theoretical distributions is the number counts of galaxies, not the mass distributions, 
the discrepancy might indicate that the simulations have too much merging in the denser 
regions. 



5. Correlation of Galaxy SED and Luminosity with Structure Location 

A number of rec ent investigations have found early-type galaxies more strongly clustered 



than the later types (ILe Fevre et al.ll2005l ; iMeneux et al.ll2006t ICoil et al.ll2006l ; ICooper et al. 
20061 ) at z = 0.5 to 2. Variation of galaxy SEDs as a function of both redshift and structure 
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location is dramatically shown in Figure [12] Here we plot the mean SED of galaxies as a 
function of z (with no selection for structures or the field) and the mean for galaxies within 
R < 1 Mpc and R = 1 - 5 Mpc from the center of each structure. This is on a scale with six 



types ranging from 1 for an E/SO galaxy SED to 6 for a starburst galaxy (IMobasher et al. 
20061 ). Figure [12] demonstrates dramatically, and in every case that the interior of the 



structures are populated with galaxies having a mean SED type lower by ~ 0.5 - 1 compared 
to the average SED type at the same redshift. iButcher fc Oemlerl ( 11984 ) first showed the 
trend for an increasing fraction of blue galaxies within clusters out to z = 0.5 and the trend 
for earlier morphol ogical types in the h ighest density regions is well known as the T - £ 
relation at z < 0.5 (IDressler et al.l 119971 ). The sample shown here for the COSMOS survey 
is the most extensive and covers a large range of redshift (0.1 < z < 1.1) using the same 
technique. An anal ysis of galaxy morp hology and environmental density in the COSMOS 
field is presented by lCapak et al.l (120061 ). Their results are consistent with those shown here. 
Figure [12] also shows a systematic gradient in the mean SED for the field galaxies - about 
+0.5 to later types from z =0.2 to 1. 

Figure [T3] shows the fields surrounding a sample of six of the LSS (#1, 2, 8, 10, 25 
and 26) with the galaxies shown in color depending on their SED type determined in the 
photometric redshift fit. (Galaxies within the Az range given in Table [3] are plotted for 
each structure.) These figures show the enhancements in galaxy density associated with 
the LSS; they also indicate the level of background contamination which any identification 
procedure must deal with. However, the most interesting feature easily seen in Figure d3] is 
the preference of the early SED-type galaxies for the denser LSS. It is important to recall 
that the sample selection used to identify the structures was involved all galaxy types, not 
just red galaxies. 



6. Evolution of Galaxy Properties with Environmental Density and Redshift 



The dependence of galaxy properties on redsh ift and environment i s one of the central 



themes of current c osmological evolution studies (|Le Fevre et al.l 120051 ; iGerke et al.l 12005 
Meneux et allbood : ICooper et allbood : ICoil et al.&od . e.g.). Here we use the overdensities 
derived as a function of the pseudo-3d space (a, 5 and z - Section [3] and Figures [3]) and 



galaxy properties (SED type, mass, luminosity and star formation rate - SFR) derived from 
broadband photometry to investigate the environmental influences. Use of the density cube 
precludes the need to identify and delineate specific LSS (Section [5]). 

Our environmental densities were derived from the surface density of all galaxies (above 
specified mass or luminosity cuts) - the densities were not derived from clusters of color 
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selected galaxies - thus the analysis below is presumably unbiased and without a priori 
correlations of environment and galaxy properties. 



6.1. Environmental Density 

To characterize the local environment of each galaxy, we use a 'relative density' measure, 
p re i defined as 

p{a, 5, z) rel = (3) 

where Ti(a,S,z) is the overdensity for each redshift slice as shown in Figures [3] and is 
the mean of this overdensity at each redshift. The mean value of the overdensity is used 
for normalization to enable comparison of widely separated redshift slices with Az = 0.1 
which have somewhat different surface densities and overdensities of galaxies due to varying 
line-of-sight depths and comoving volumes. (The relative densities p may be translated back 
to E (Mpc- 2 per 0.1 in z) using E(z) = 1.2, 0.52, 0.29, 0.17, 0.14, 0.19, 0.18, 0.15, 0.13 
and 0.06 for z = 0.15 to 1.05, sampled every 0.1 in z.) For each galaxy (Section 12.11) . the 
environmental density was obtained from the pseudo 3-d cube using its a, 5 and best fit 
photometric redshift. 



6.2. Sample Selection and Completeness 

It is of course vital that the sample selection function (see Sections 12.11 & 12.21) not intro- 
duce biases as a function of redshift which masquerade as changes in the galaxy properties. 
We make use of two alternative samples with : #1) a mass cutoff of M* > 3 x 10 9 M (82,274 
galaxies) and alternatively, #2) a luminosity cutoff with My < —19 mag (101,018 galaxies). 
Figures [2] show the observed distributions of My and M* as a function of z. As discussed in 
Section I2.2[ there is little change in the mass function at z < 1 and therefore most of the 
variations in the mass function at M < 5 x 10 9 Mp , are likely the result of incompleteness 



at z > 0.7 (see Figure [2J right panel). iBorch et al.1 (120061 ) found a possible doubling of the 



integrated mass function of galaxies from z = 1 down to 0.2 (in the COMBO-17 survey). We 
take this as an upper limit since the sample used here shows no significant variation aside 
from the aforementioned incompleteness (see Figure [2]- right). Similarly, the selection on My 
is chosen to be close to the limit at which completeness starts to become an issue (see Figure 
Uleft). 



We develop the two samples in parallel since one cannot assume a priori that the galactic 
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masses and/or luminosities are invariant from z = 1 to 0. For example, one expects the 
luminosity of each galaxy to vary at z <1 (even in the absence of further star formation 
or merging) due to dimming as the stellar population ages. For this reason, adoption of a 
fixed My cut would yield a sample with larger surface density at z =1 than at z =0.2. The 
fixed mass-cut sample most likely comes closest to generating equivalent galaxy samples at 
z < 1.1; however, at higher redshifts, it is likely the masses will be changing more rapidly. 
A later paper will explore various evolution scenarios for the luminosity-sleeted sample. 

More conservative higher mass and luminosity cutoffs would of course yield greater 
completeness at high z. On the other hand, since the various basic galaxy types have quite 
different masses and luminosities, this would compromise one's ability to probe evolution 
between types. Specifically, a very high mass cutoff (e.g. M >fewxl0 10 M ) would largely 
limit the samples to just the most massive E's and spirals and under-represent the lower 
mass, late type systems which have significant star formation activity. This would severely 
compromise the dynamic range that could be investigated vis-a-vis the transformation from 
late to early type galaxies. 

The distribution of p for the sample of galaxies was also examined to select a lower 
cutoff in density for the analysis. This was required since very low overdensities, compared 
to the mean background, are not quantitatively meaningful - in areas where the adaptive 
smoothing detects no significant overdensity exceeding 3a (see Appendix |A|) , it smoothes 
the surface density down to a value detemined by the largest spatial-smoothing width. The 
adopted density cutoff reduced the final samples to 10,382 and 12,523 galaxies for samples 
#1 and 2, respectively. (The lowest overdensity to which one may carry this analysis is 
determined by the background counts of galaxies at each redshift. This is, in turn, largely 
a function of the photometric redshift accuracy. Higher accuracy photometric redshifts will 
enable extension of this investigation to lower density and the field.) 



6.3. Galaxy Properties : Mass (M*), SED type, Early- Type Fraction, My, 

SFR and t sf 



Galaxy SED types and rest-frame luminosities (My) are by-products of the photometric 
redshift fitting. Their masses were derived using the intrinsic SED to estimate the mass- 
to-light ratio together with the absolute V magnitude obtained from the observed fluxes 
(IMobasher et al.ll2006l ). The SED types range from 1 to 6 with : 1 = E, 2 = Sa/Sb, 3 = 
Sc, 4 = Im, 5,6 = two starburst populations (defined by lKinney et al.lll996l ). The early-type 
galaxy fraction was calculated, taking all with SED type < 1.9 to be 'early-type'. For each 
galaxy, the SFR was estimated from the intrinsic SED and observed fluxes, extrapolated 



into the UV. (As with the mass estimates, the SFRs have been aperture-corrected using 
the auto-magnitude parameter from SEXTRACTOR.) We use the SFR estimated from the 
extinction-corrected, rest-frame A =1500A continuum ( Mobasherl 120061 ) . 



We also calculate the ratio of the galaxy mass to the SFR, yielding a characteristic 
timescale to form the existing galactic mass of stars at the currently observed SFR (specifi- 
cally, tsf = M*/ SFR). For a starburst t$f will be significantly less than the Hubble time at 
the observed redshift whereas a galaxy for which the current star formation is relatively low, 
compared to that in the past, will have a long t sf . t S f is equal to the inverse of the specific 
SFR per unit stellar mass of the galaxy, sometimes called the 'star formation efficiency' 
(SFE). 

The galaxy samples were binned using 4 equal z bins of width Az = 0.23 between z = 
0.2 and 1.1 and 4 logarithmically spaced bins in density p from 8 to 215. For the adopted 
cosmology, the redshift bins are centered at lookback times of ~ 3.5, 5.3, 6.6, and 7.7 Gyr. 
Within each bin, the median values of each galactic property were determined. The median 
was used rather than the mean since it is less susceptible to a few extreme values and hence 
the uncertainty in the median estimates can be small even for samples with a large intrinsic 
dispersion. To estimate the uncertainties in the median values, Monte Carlo simulations were 
done on the observed distributions, adding randomly sampled uncertainties from a normal 
distribution. We adopted uncertainties (la) in each of the bolometric quantities (M*, My 
and SFR) of a factor of 2 from their nominal values for each galaxy; for the SED type, we 
assume an uncertainty of ±1 for the type. (The factor of 2 uncertainty is an approximation 
to allow for uncertainties in photometric calibrations and the SED fitting.) The median was 
measured for each of 500 simulations, and the dispersion of the median distribution was 
taken as the uncertainty in the median for the observed sample. 



6.4. Galaxy Evolution with Redshift and Density 

In Figures [T41 to [T§1 the median stellar mass M#, My, SED, SFR and t sf are plotted 
for each redshift range as a function of density. The SEDs, SFR and t$f all exhibit very 
significant variation as function of both redshift and density. The mass and luminosity 
distributions are partially affected by selection bias but only in the highest z bin. 

The median mass and My distributions (Figures dH and [To]) and comparison of the 
mass- and luminosity- limited samples (left and right panels) may be used to assess the 
possible influence of incompleteness at the highest redshifts. The median masses show no 
systematic increase with z except in the z = 0.88 to 1.1 bin for which the masses appear 
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systematically higher by a factor of 1.5 to 2 compared to lower z. This increase is very 
likely due to sample incompleteness or Malmquist bias since the sample is deficient in the 
galaxies with M # < 5 x 10 9 M (see Figure [2], right panel). My exhibits a somewhat larger 
(AMy ~ 1 mag) increase - some of this is probably also due to incompleteness (see Figure 
[2] left panel), but since it is larger than the mass shift, some of the My variation is likely due 
to actual evolution of My in the g alaxies. Passive evo lution of the stellar populations from 



z = 1 to is AMy ~ 1.2 mag (e.g. iDahlen et al.ll2005l ). To summarize, modest variations in 



the mass and luminosity medians between z = 0.8 and 1 are probably due to incompleteness; 
at lower redshifts, no such variations are seen and the samples are probably complete to 
better than ~ 20%. 

The median masses clearly grow with increasing density, at each redshift. This cannot 
be a sample selection bias since that should be constant at each redshift. The doubling 
of the median mass in high density environments compared to lower densities is seen at 
all redshifts out to z = 1; this is undoubtedly reflecting the fact that the early type SEDs 
also are more prevalent in the denser environments (see below) and these are often massive 
galaxies. Consistent with this interpretation is the fact that although M* (Figure fT4"l) exhibits 
dependence on p, My does not (Figure [T5]), implying that the median mass-to-light ratio is 
lower at high density. 



6.4-1- Galaxy Spectral Type & Early-Type Fraction 



The galaxy SED types shown in Figure [16] exhibit very significant variations with both 
redshift and density - in the sense that earlier type SEDs (E's) are seen at higher density, 
later types in the lower density regions. And for all densities, the median galaxy type is later 
(i.e. bluer, star forming) at higher redshifts. The major variation with z occurs between the 
lowest two redshifts (from z ~ 0.3 to 0.5 or lookback times less than 5.3 Gyr) - all three 
high z bins have similar SEDs and their variations with density are the same. Numerous 
studies have noted the strong increase in the fraction of early type galaxies (classified by 



both SEDs and morpho l ogies) in dense envir o ments out to z > 1 (e.g. iDavis fc Geller 



1976 



Postman k Gelled 11984 ; iDressler et all Il997l : iKodama et all 12004 ; iKauffmann et all 120041 ) . 
The results shown in Figure [16] show very clearly that similar density correlations are seen 
over the entire redshift range. The mean SED shifts to earlier type at lower z for both low 
and high density environments. The actual surface density for the breakpoint between the 
late and early SEDs shifts does not appear to shift more than a factor of 2 since the break 
occurs at approximately the same relative density p and the density normalization from z 
~ 0.3 to the higher z's changes by less than a factor of 2 (see Section [6TTT) . 
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The percentage of galaxies with SED type < 1.9 is shown in Figure [T71 exhibiting 
variations like those in the median SED type (as it should, since they are closely related). 
We include the early type fraction since it is often used to characterize the galaxy populations 
in evolutionary studies. As with the SEDs, the major shift with z occurs beween the lowest 
two redshift bins, and at all redshifts an increased early-type fraction in seen above p ~ 50. 
Figure [17] clearly demonstrates that the galaxy- type correlation with density was clearly in 
place before z = 1 and we have extended this correlation to low densities, as well as the 
dense clusters. 



6.4-2. Star Formation Rates and Timescales 



The median SFRs per galaxy (Figure [T8j) rise systematically with redshift for all densi- 
ties. The overall increas e by a factor of 4.5 f rom z = 0.3 to 1 is similar to that found in many 



earlier studies (e.g. see iMadau et al.lll996l; iHopkind |2004 lJuneau et al.l 120051 ; iBundy et al. 



20051 ; iBell et al.ll2005l ; ISchiminovich et al.ll2005l ). The observed increase at higher redshift is 
extremely well fit by a linear dependence on lookback time over this range, Ti 00 kback = 3.5 
to 7.7 Gyr. Figure [T8l also shows evidence of a slight decrease in the median SFR at higher 
densities, with this decrease being steepest at low redshift (z~ 0.2). The steep decline in 
the SFR to lower redshift is possibly due to the depletion of ISM to fuel star formation 
and AGN/SF feedback processes. Discriminating between these may be accomplished with 
future observations of the star forming gas content with the ALMA array. 

Normalizing the SFRs by the stellar mass of each galaxy, the SF timescale (tsf, Figure 
[T9~l) shows much stronger density correlation than the SFR. At all densities, the SF timescale 
is a factor of 2 - 3 shorter in all three high redshift bins compared with z = 0.2 to 0.43. 
And a factor of 4 - 5 increase in the SF timescale occurs between the low and high density 
environments at all redshifts with the strongest density dependence occurring at the lowest 
redshift. These results imply that most of the stellar mass in dense environments must have 
formed much earlier than z = 1 whereas a significant amount (~25%) of the stellar mass in 
the low density environments must have formed at z = 1.1 to 0.4 (based on the measured 
SF timescales). 



6-4-3. Downsizing of Star Forming Galaxies - the Maturity Parameter (fi) 



A number of investigations have suggested that star formation occurs earlier in the 
most massive galaxies and as the universe ages the star formation progresses to less and 



-24- 



less massive systems , a phenomenon of ten referred to as 'downsizing' (ICowie et al.l 11996 



Kodama et al.l I2004J ; iBundy et al.ll2005l ). However, this phenomenon can be blurred and 



sometimes confused with the earlier formation times for galaxies in dense clusters coupled 
with the high abundance of massive galaxies in clusters. Here, we attempt to separate these 
affects to investigate the relative formation times of high versus moderate mass galaxies as 
a function of both redshift and density. 

For this discussion, we define a parameter which we will call the Maturity (/i), equal to 
the ratio of the star formation timescale (tsf used above) to the cosmic time (r cosrnic = the 
age of the universe at each galaxy's redshift). With this definition, the Maturity is unity if 
the observed stellar mass could have formed at the observed star formation rate within the 
age of the universe (at the redshift of the galaxy). The Maturity will be < 1 (youth) if it 
is forming stars at a sufficiently high rate that its mass could be produced in less than the 
cosmic time; the Maturity will be > 1 (middle to old age) if its current SF rate is low and 
most of its stars must have been formed earlier (with \i <1, youth) at a star formation rate 
much higher than that presently measured. Obviously, initial starburst systems would have \x 
<< 1 and old elliptical galaxies \i » 1. The Maturity, defined in this manner, will continue 
to increase at later cosmic epochs if the star formation rate remains low. On the other hand, 
if the aging galaxy undergoes a late-life starburst (mid-life crisis), it will be rejuvenated (/i 
Jj). But, if the the starburst is brief and not substantial, the galaxy will return more or less 
to its prior state of Maturity after the starburst. (As with humans, rejuvenation may be 
superficial and illusory ! Our use of medians for charting the overall evolution of the galaxy 
populations probes the typical Maturity, thus avoiding the 'noise' due to short starbursts. 
[Anthropomorphizing this galactic parameter can actually help to visualize and track the 
galactic changes associated with evolution of /i.] (T his Matur i ty par ameter is similar but 
not identical to the 'Birthrate' quantity discussed by iBell et al.l (120051 ). but 'maturity' more 
aptly connotes what this parameter characterizes. ) 

In Figure [20] the Maturity is shown as a function of both redshift and density, separately 
for galaxies of high and low mass. The two samp les were separated at M* = 5 x 10 10 M Q . 
Kauffmann et al.l (120041 ) and iKodama et al.l (120041 ) found a 'break' between old, red galaxies 
and younger, blue galaxies at a mass of ~ 3 x 10 10 M o at z = and 1, respectively. We 
have adopted a s omewhat higher valu e in order to minimize incompleteness in the low mass 
sample at z > 1. IBundy et al.l (120051 ) argue that the break mass varies with redshift, rising 
to ~ 10 11 M G at z =1; however, we do not see such clear variation (see below). In any case, 
we have found by experimenting with the mass cut that a factor of 2 variation in the vale 
mass cut (from 5 x 10 10 M & ) did not change the behaviors discussed below. 



Figure 1201 shows that at all redshifts and densities probed here, the more massive galaxies 
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are always more mature than the lower mass galaxies. At each redshift and environmental 
density, the lower mass galaxies are systematically 5-10 times less 'mature' than the massive 
galaxies. Once again, we emphasize that these are not color-differentiated galaxy samples 
- just mass-differentiated for which there is no a priori association with 'age'. Although, if 
the more massive galaxies tend to be more mature, obviously, they will appear redder. 

For both the high and low mass galaxies, the median Maturity is either constant or 
increases with time, i.e. to lower redshift (Figure [2171) ; it never decreases with time - as 
well it could if the star formation in either environment was delayed to commence at a 
late epoch (e.g. some dwarf galaxies). A constant Maturity implies a steady SFR over 
cosmic time, whereas an increasing maturity suggests diminishing SF with time. The more 
mas sive galaxies clearly must have had an early phase of rapid star formation at z < 2 



(c.f. lJuneau et al.l 120051 ) with relatively little star formation at z < 1.1 in order to appear 
so mature (// ~ 1 to 2 at z ~ 1, Figure [201) . By contrast, the lower mass galaxies exhibit 
fairly constant immaturity down to z ~ 0.43, implying that on-going, fairly constant star 
formation has occurred from z = 1.1 to 0.43 and very likely also at the high z. However, at 
the z < 0.43 the maturity of the lower mass galaxies rises in all environments, implying a 
significantly decreased star formation rate for lookback times < 5 Gyr. 

Figure [20] suggests that galaxies of all masses (at z < 1.1) are more mature in the 
dense environments, not just the high mass galaxies! The lowest redshift bins both show 
T s f/r cosmic rising quite significantly at the highest density while at the other redshifts, a 
factor ~ 2 increase in the maturity occurs between the lowest to highest densities. This 
clearly requires that the epoch of rapid star formation for galaxies of both high and lower 
mass must be earlier in the denser environments than in the field. 

In Figure [5T] the Maturity is shown as a function of galactic mass for high and low 
density environments, with separate plots for each redshift. The cut between high and low 
density was taken at p — 45, i.e. between the middle two bins of the 4 density bins used 
earlier. In Figure [21], the overall distribution of galaxies for all densities is shown by the 
colored shading while the high and low density environments are shown in the red and blue 
contours, respectively. The separation of the old (/i> 1) and young ( /i< 1) galaxies is 
seen as a bimodal distribution and their loci change systematically with redshift. Evidence 
of evolution i n the mass separating starfor m ing and non-s t ar form ing has been claimed by 



Bundv et al.l (120051 ) and iBorch et al.l (120061 ). iBundv et all (120051 ) found M break oc (1 + z 



V'l. 



however, we find it difficult to identify a distinct mass which can be said to divide the 
mature and immature populations, since at most masses between 10 10 and 10 11 M , /i can 
range from < 1 to > 10 (see Figure |2~T]) . Most of the mature galaxies occur above fewxlO 10 
M Q , but there is no sharp cutoff at most redshifts. In fact, for the lowest redshift bin, two 
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distinct, mature sequences can be seen at 5 and 20 - 30. The latter could correspond to 
the maturation of the /i~ 6 to 8 sequence seen at higher z; the former might correspond to 
maturation of the fi< 1 galaxies seen at earlier epochs. Possibly, the combination of these 
two mature seq uences at low redshifts account for the apparent evolution of the break mass 
as discussed by iBundy et al.l (120051 ). Future work is planned using the COSMOS GALEX 
UV measurement to verify this second mature sequence. 

Galactic down-sizing with the most massive galaxies forming earliest is of course at 
variance with the expe ctation of the most simpl e hiera rchical galaxy formation scenarios. 
However, iKaiserl (119841 ) and later ICen fc Ostrikerl (119931 ) suggested a model of biased galaxy 
formation with the most massive galaxies forming within the highest peaks of the initial 
density field, and this is a commonly accepted explanation. In the highest peaks, there is 
more mass available for buildup of the most massive galaxies and th e rate of growth is higher 
where the density of sub-halos is higher (e.g. iDe Lucia et al.ll2004l ). The results shown here 
suggest that even at relatively low environmental densities, the more massive galaxies are 
formed earlier than the low mass galaxies - although not as early as the massive galaxies 
in very dense environments. This suggests that the formation of massive galaxies occurs by 
two processes - one local, responsible for the early growth of massive galaxies in low density 
regions, the other associated with high overdensity regions where the growth occurs more 
rapidly and in some cases is carried to the very highest galactic masses. 



7. Summary and Conclusions 

The COSMOS photometric redshifts now have sufficient accuracy [a z j(\ + z) — 0.03) 
to enable identification of LSS at z = 0.1 to 1.1. We have developed an adaptive smoothing 
procedure to be applied to the galaxy density distributions in photometric redshift slices with 
Az = 0.1 to identify LSS on scales less than 1 Mpc up to 30 Mpc with optimal signal-to- 
noise ratio across the range of spatial scales. This procedure has been tested with excellent 
results on mock redshift slices and on the dark matter particle distribution from a ACDM 
simulation (see Appendices IA.2I and IA.3j) . 

The adaptive smoothing is applied to the COSMOS photometric redshift catalog with 
selection z < 1.1, 19 <Iab <25 mag and M v < —18 mag - a sample of 150,000 galaxies. No 
color or SED selection is imposed, so that the defined structures are intrinsically unbiased 
with respect to galaxy type. From the galaxy over-densities derived from the adaptive 
smoothing, we have delineated 42 LSS and galaxy clusters in the pseudo 3-d space (a,S,z). 
The surface density plots of the structures are shown in Figures H] and [5] their measured 
properties are given in Table [3J Five of the most massive structures have stellar masses 
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(determined from the galaxy photometry) of M* > 10 M Q . Several have extents which can 
be traced over 10 Mpc (comoving). Their total masses including dark matter are likely to be 
50 - 100 times greater. The Richness of the core regions of these structures is typical of Abell 
class 1-3. The derived mass function for the LSS is consistent with the total mass function 



for clusters derived by (jBahcall fc Cenl Il993l ; iReiprich fc Bohringerl |2002| ) from optical and 
X-ray studies. 



a ( 


2006) and 


Cassata 


f 


2006) 



emissi on, many of which are located within the LSS discussed here, are discussed by lFinoguenov et al. 



(120061 ). These clusters are identified optically by wavelet analysis of the early type galaxies 
in the COSMOS photometric redshift catalog and from the diffuse X-ray emission. We have 
compared the fractional areas seen at different overdensities and find general agreement to 
within ~ 50% with the predictions of ACDM simulations (processed similarly) - with less 
than 1% of the areas of the redshift slices having overdensities exceeding 10:1. However, the 
observed filling factor distribution does reach higher overdensity and this may indicate that 
the simulations have too high an efficiency for merging in dense regions., 

We have investigated the dependence of galaxy evolution on environment using the 
structures defined here and the SED types taken from the photometric redshift fitting. We 
find that in every structure the mean galaxy SED type within the high density core of 
the structures is earlier than the mean SED type at the same redshift. Our study thus 
confirms, with a sample of 4 2 structures/clusters , the correlation of galaxy evolution with 



environmenta l loca tion (e.g. iDressler et al.l 119971 : ISmith et al.l 120051 : iPostman et al.l 12005 



Cooper et al. 


2006 


to 1. 


Capak et al. 



and r eferences cited therein)) over the full range of redshift z = 0.1 
(120061 ) find a similar result using an entirely independent measure of 



environmental density and using galaxy morphology instead of SED type. 

Extensive analysis was done to analyze the correlations of galaxy properties (SED, mass, 
luminosity and SFR) with redshift and enviroment. The median SED type and star formation 
activity varies strongly with both redshift and environmental density. The maturity of the 
stellar populations and the 'downsizing' of SF in galaxies are both strongly varying with 
epoch and environment. Although the more massive galaxies clearly tend to have lower 
SFR per unit galactic mass, we question whether it is possible to define a distinct 'break 
mass' separating active and inactive star-forming galaxy populations. And over the range z 
< 1.1, we don't see strong evidence of evolution in the masses of galaxies undergoing active 
star- formation (at the level of < a factor 2). 
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A. LSS Identification with Adaptive Smoothing 



In the past, a number of algorithms or techniques have been used for automated identi- 
fication and cha racterization of galaxy c lustering, including percolatio n and Voronoi tessela- 
tion techniques (van de Weygaert 1994; lEbeling fc Wiedenmann 1993 ; Marinoni et al. 120021; 



Gerke et al.ll2005l ). wavelet analysis (lEsca 



and matched filter (IPostman et al.lll996 



era &; M acGilli vrav 



1995 



Schuecker &; Boehringer 



Finoguenov et al.ll200q ) 



19981 ). An algorithm for 



the identification of structures must be capable of detecting structures on multiple angular 
scales, and with only low order assumptions regarding the internal density profile of the 
structures. Techniques which search for a particular scale or assume, a priori, a density 
profile (or equivalently a spatial weighting function) will have highest sensitivity for struc- 
tures with the specified parameters, thereby biasing a derived distribution function for the 
recovered structures. It is also highly desirable that the algorithm be capable of displaying 
compact structures simultaneously with more extended, low density structures. Presumably, 
within large structures there will be high density substructures which one would not want 
to smooth out into low spatial frequencies. Conversely, if a high density structure is fully 
detected at high spatial frequencies, one would not want its power to be carried out to low 
spatial frequencies as an extended halo. Multi-scale algorithms like wavelet and adaptive 
smoothing seem therefore most appropriate. 



For the structure identification, we have developed an adaptive kernel smoothing algo- 
rithm, specifically tailored to have these characteristics. 
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A.l. Algorithm 

The algorithm consists of a loop, starting at low smoothing width, going to successively 
larger smoothing kernels, removing power from the current 2-d residual 'image' if it exceeds a 
specified signal-to-noise ratio at the current level of smoothing. The 'image' being processed 
is the projected surface-density (E) of galaxies in a redshift slice. Starting at the initial 
highest resolution (n=l), we calculate the smoothed surface- density (E n ) and background 
surface-density (B) from 

E n = E' n _ x * K n (Ala) 

B n = (E;_ x - A n ) * K 2n (Alb) 

where V is the convolution operator, K n is a 2-d smoothing kernel of width n, normalized 
such that its integral is unity and the Kernel K 2n used to convolve the background has twice 
the width, ie. 2n. The 'power available' at resolution n is calculated as 

A n = E n - B n (A2) 

Since Eq. lAlbl depends on IA21 these equations are iterated (typically 4 times) to arrive at 
the 'best' estimates of the background (without the high frequency power included) and the 
A„ with the most low frequency, background removed. 

If a n is the noise 'image' at resolution n, then the signal-to-noise ratio, S(A) n , on the 
delta residual-density is then 

S(A) n = ^ (A3) 

On 

and the signal-to-noise ratio, S(E)„, on the original surface- density, smoothed to resolution 
n, is 

5(S) n = (A4) 

On 

The power to be removed at resolution n is then given by 

P n = {H {3{A) n - snr 5 ) ■ H {S{T) n - snr total )} ■ A n (A5) 

where H(x) is the Heavyside function (H=0 for x < 0, H=l for x > 0). 'snr^' is an adjustable 
parameter specifying the minimum signal-to-noise ratio in the residual image required before 
power is removed at width 'n'. Similarly snr to4a ; is a parameter specifying the minimum 
signal-to-noise ratio required when the original total-power surface-density is smoothed to 
resolution 'n'. Having these two conditions is crucial to the excellent results obtained with 
this procedure - allowing small values of snr^ to be used while avoiding the retention of 'noise' 
peaks. For any pixels which do not satisfy this double criteria for signal-to-noise ratio, the 
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residual power is retained to the next level of smoothing. The residual image (with lower 
spatial frequency power) to be used as input on the next iteration at larger smoothing kernel 
(n+1) is therefore given by 

K = K-l - Pn (A6) 

Steps Al to A4 are repeated with successively larger values of V up to n max . 

After reaching ia max , the adaptive smoothed surface-density (X/maz) is then given by 

V final = K + Y, P » ( A? ) 

n 

The procedure described above has the following desirable features : 

1) It is conservative, i.e. the 2-d integral of the original and final surface- densities are equal. 

2) Power is retained at the highest spatial frequencies and not smoothed out to lower fre- 
quency as long as its signal-to-noise ratio is sufficient (i.e. greater than the specified snr^). 

3) High frequency power is removed first, extended haloing around high density regions is 
thus minimized. 

4) Features seen in the final adaptively smoothed surface-density have a well-determined 
significance and resolution. 

One caution : since the resolution is variable across the adaptively smoothed 'image', 
the usual intuition that judges significance or signal-to- noise ratio by comparison with the 
amplitude of high-frequency noise is not reliable. 

There are several parameters which are important to results of the adaptive smoothing 
process outlined above : 

1) the signal-to-noise ratio, snr tota i, to be required in the total surface- density, S (smoothed 
to the current resolution). This parameter is set at snr tota ; = 3 so that virtually all features 
seen in the final adaptively smoothed image will be 'statistically' significant. 

2) the signal-to-noise ratio, snr,5, specifying whether the A n signal is removed before 
proceding to a lower resolution filter. This parameter should be set such that power is 
removed at the highest spatial frequency for which there is a 'reliable' signal, but avoiding 
removal of what is essentially low-frequency power before 'its time has come'. Based on trial 
and error, we have adopted snr,5 =1 for the LSS identification. Although it might seem that 
1 a would be risky, the 3 a condition (above) assures that most features will be significant. 

3) the maximum filter width, n max . The maximum filter width was taken at 0.33°, i.e. 
23% of the linear size of the COSMOS field. 
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Two smoothing kernels were used : a boxcar and Gaussian. The boxcar was used 
for program development since it was faster, but all final results employ a 2-d, symmetric 
Gaussian filter (implemented with a Fourier transform for speed). It is well known that 
boxcar filters can introduce high spatial frequency edges whereas the Gaussian is better 
behaved in this respect. 

A. 2. Simulation Tests 

To test the algorithm described above with conditions similar to the galaxy counts in 
the COSMOS photometric redshift catalog, we have simulated a single redshift slice with 
10,000 galaxies. Approximately half of the galaxies were distributed within Gaussian profile 
structures with a distribution of peak densities and sizes. The other half of the sample 
galaxies were distributed randomly across the field. 

The parameters for the Gaussian-profile structures in the simulation are listed in Table 
[5j Figure l22h shows the input surface-densities profiles. For the lowest 3 rows, the simulated 
structures have increasing peak surface-densities toward the top and increasing in size going 
to the left. The top row simulates more complex structures with three internal components 
having varying sizes and surface densities. Figure [22b shows the simulated distribution 
of galaxies, consisting of 5260 galaxies randomly placed and 4740 galaxies populated with 
probability given by Figure [22k . It is important to realize that since the density profiles 
of the structures are sampled randomly, the simulation distribution, input to the adaptive 
smoothing algorithm, is not identical to that shown in Figure [22b . i.e. there is shot noise. 
Therefore, the algorithm should not be expected to return the smooth input distributions 
(Figure [22k) exactly. 

Figure [23] shows the surface- density recovered from the galaxy distribution shown in 
Figure [22b using the adaptive smoothing. In fact, the algorithm has done an excellent job of 
recovering all structures which were statistically significant in Figure [22b . including the top 
row with complex, internal structure. The three structures in the lower right of Figure [22k 
were not recovered but these were all sufficiently low in surface- density and/or size that their 
total numbers of galaxies were not statistically significant (2, 8 and 6 galaxies respectively - 
see Table [5]) . Lastly, it is worthwhile to emphasize that the algorithm did not find structures 
which were not in the input simulation, i.e. noise in the random galaxy population was not 
falsely detected using parameters for the simulation distribution (numbers of galaxies and 
fraction in structures) and for the detection algorithm similar to those used for the COSMOS 
structure detection. 
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A.3. Test on ACDM Simulations 



As an additional test we have ap plied the adaptive s moothing; procedure to one of the 



Virgo consortium ACDM si mulations (IBenson et al. 



simulation COSMOS wedge (jSpringel et al 



2005 



2001) and the more recent Millennium 



Croton et al.l 120061 ) 



The Virgo simulation had dark matter particles of mass 1.4 x 10 n /h M Q and for the 
purposes of our test, we sampled the dark matter particles to obtain a surface density of 
particles in each redshift slice similar to that of galaxies in the COSMOS photo-z catalog. 
This was done to keep the simulation shot noise characteristics similar to those of the obser- 
vational data being analysed here. The results for adaptive smoothing of the ACDM Virgo 
simulation are shown for z = 0.35 and 0.93 in Figure [2H The algorithm reliably recovers 
all significant structures seen in the simulation. It is noteworthy also that the scale of the 
structures seen here is qualitatively similar to that actually found in our application to the 
COSMOS photo-z catalog. Compare Figure EH with the similar redshift frames of Figure HI 

In our tests on the Millennium Simulation, the objective was to determine if similar 
structures were seen in this most up to date simulation as in the COSMOS data. Thus 
the galaxies in the Millennium COSMOS light cone were each given a redshift uncertainty 
similar to that in the COSMOS photometric redshifts and then processed identically to the 
observed galaxies. In Figure [25] the overdensities from the simulation are integrated along 
the line of sight from z = 0.25 to 1.05 as was done for Figure including keeping the 
same contours for both plots. Extremely good correspondence is seen from the comparison 
indicating that the adaptive smoothing is recovering very similar structures in both and by 
implication, both theory and observations have similar intrinsic structure. In Section HI we 
make a more quantitative comparison by measuring the area filling factor as a function of 
overdensity and redshift. 
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Table 1. LSS Galaxy Samples 



Sample 


z I cutoff (mag) 


# of galaxies 


# (early type) a 


# (late type) a 


all-z 


z < 3. 26 


487973 


59285 


428688 




25 


228487 


34303 


194184 




24 


106487 


25278 


81209 


low-z 


z < 1.1 25 


150162 


25583 


124579 



Note. — All samples have M.y < —18 and galaxies must be detected in at least 4 bands. 
a Early type : SED type < 2.5; Late type : SED type > 2.5 



Table 2. Galaxy Selection 'Completeness' 



Redshift interval 


Rel. # per unit vol. 


Rel. Mass per unit vol. 


0.1 - 0.3 


0.52 


0.69 


0.3 - 0.5 


0.69 


0.94 


0.5 - 0.7 


1.0 


0.92 


0.7 - 0.9 


0.56 


0.96 


0.9 - 1.1 


0.47 


1.0 



Note. — The total number of galaxies and total mass of galaxies 
per unit comoving volume were evaluated by integrating the distribution 
functions shown in right panel of Figure [2] at masses above 10 9 M©. 
These totals were normalized to the values in the redshift interval with 
the largest values in order to assess the count and mass incompleteness 
relative to this maximum bin. 
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Table 3: Structures in the COSMOS field 



FWHM b 



Structure # 

f"X" rav id'!* 1 


RA-150 a 


Dec-2 a 




ARA 


ADEC 


Az 
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ze c 


% LSS d 


galaxies 


M,' 


Central ^XAlpc 


r°i 

I 1 


\ ) 




(°) 


(°) 






Mpc 






10 rz M Q 




± y t t , ±VU, ±Uo, ±UO ) 












22 





17 





27 


12.69 


35 


1767 


23.63 


188 




15 


20 


88 





20 





21 





25 


17.41 


30 


8 1 5 


15.72 


'■>'■) 




-0 33 
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18 





17 





2(1 
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■31) 
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4S 




0.57 
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0.20 
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00 
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01 
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0.20 





01 
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00 
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70 







01 




39 


0.33 


-0.36 


0.60 





01 





01 





00 





63 


36 







01 


10 


40 
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0.31 


0.30 
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-0.30 


0.51 
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00 





38 
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a RA and DEC (J2000) of peak galaxy surface-density and the centroid z. 

b Full- width at half-maximum evaluated from the 2.3 xc where a is the dispersion (from the calculated 2nd 
moment). 

Estimated as (ARA 2 + ADEC 2 ) 1 / 2 converted to co-moving Mpc. 

d Mean probability that a galaxy within the structure is within the structure rather than being in the projected 
background, estimated as Y^lss/^lss + ^Background)- 

e Total number of galaxies estimated within structure, corrected for 'field or background' contamination by 
comparing the surface-density in the structure with the background surface-density for that redshift slice - 
see note (d) 

•^The total stellar mass (M„), estimated from the absolute magnitude of each galaxy and using a mass-to-light 
ratio appropriate to the galaxy SED. For each galaxy, this photometric mass is multiplied by the probability 
that it is within the structure rather than being a projected field/background galaxy (see note (e)). 

s The central surface density used to evaluate the structure/cluster Richness - the number of galaxies within 
radius 1 Mpc brighter than 2 magnitudes below the 3rd brightest galaxy. This column is blank if there are 
too few galaxies to estimate a Richness (i.e < 10) 

"In parenthesis we give the cross-reference to the ID for the X-ray clusters from lFinoguenov et al.1 ([2006). The 
wavelet technique used bv lFinoguenov et al.l (|2006T ) on the X-ray emission and on early SED type galaxies in 
the photometric redshift catalog is selective toward more compact structures than the adaptive smoothing 
technique used here; therefore, in many cases, these cross-identifications should be viewed only as 'possible', 
based on close proximity in a, S, and z. 
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Table 4. Relative Variances for 2.5 deg 2 field 



K cdshift 


Comoving Volume 


M(vol) 


cr(M(vol)) 


Halo mass 


Number 


< b > 


"shot 


^cosmic 


a(M,z) 


range 


10 6 present-day Mpc 3 


10 17 M e 




range (M Q ) 


expected 










0.2 - 0.4 


0.80 


0.327 


0.149 


10 13 -10 14 


142 


1.62 


0.084 


0.240 


0.254 


0.2 - 0.4 


0.80 


0.327 


0.149 


10 14 -10 15 


6 


3.03 


0.397 


0.451 


0.601 


0.4 - 0.6 


1.78 


0.727 


0.099 


10 13 -10 14 


279 


1.84 


0.060 


0.182 


0.192 


0.4 - 0.6 


1.78 


0.727 


0.099 


10 14 -10 15 


10 


3.52 


0.324 


0.349 


0.476 


0.6 - 0.8 


2.79 


1.14 


0.075 


10 13 -10 14 


372 


2.09 


0.052 


0.157 


0.165 


0.6 - 0.8 


2.79 


1.14 


0.075 


10 14 -10 15 


9 


4.08 


0.327 


0.307 


0.448 


0.8 - 1 


3.67 


1.50 


0.061 


10 13 -10 14 


405 


2.37 


0.050 


0.145 


0.153 


0.8 - 1 


3.67 


1.50 


0.061 


10 14 -10 15 


7 


4.71 


0.371 


0.288 


0.470 


1 - 1.2 


4.40 


1.79 


0.052 


10 13 -10 14 


387 


2.69 


0.051 


0.140 


0.149 


1 - 1.2 


4.40 


1.79 


0.052 


10 i4_ 10 i5 


5 


5.41 


0.458 


0.282 


0.538 



Note. — The field size of 2.5 deg 2 was taken to match the area of the sample taken from the COSMOS photometric rcdshift catalog 
used here. The table includes the expected number of halos in two mass ranges, together with the relative variances due to shot noise, 
cosmic variance and the combined total variance. Values calculated for cosmological parameters : Ho = 70 km s _1 Mpc -1 , Um = 0.3, 
fl A = 0.7 and ct 8 = 0.74. 
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Table 5. Parameters for Structures in Simulation 



Structure # 


RA 


Dec 


FWHM(') a 


Peak (#/D°) b 


# galaxies 


1 


-0.48 


-0.64 


1.0 


6148 


2 


2 


-0.16 


-0.64 


1.9 


6148 


8 


3 


0.16 


-0.64 


5.8 


6148 


65 


4 


0.48 


-0.64 


11.5 


6148 


252 


5 


-0.48 


-0.24 


1.0 


18446 


6 


6 


-0.16 


-0.24 
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Note. — Total number of galaxies in all the stuctures is 4740. 
a Full-width at half maximum (arcmin) for the Gaussian galaxy distribution. 
b Peak surface density of galaxies. 
c Total number of galaxies in structure. 
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Fig. 1. — a) Lef t : The number distr ibution of galaxies from the COSMOS photometric 
redshift catalog (IMobasher et al.ll2006l ) is shown as a function of redshift in bins of width 
Az = 0.2 for different I-band magnitude cutoffs, b) Right : The widths of the derived 
photometric redshift probability distributions are shown as a function of z. The quantity 
plotted is the mean of 1.3xa z / (1 + z) for the photo-z likelihood distributions for all galaxies 
in each redshift bin of Az = 0.2. For a Gaussian uncertainty distribution, 50% of the 
probability is within a full width of 1.3 xa z of the peak. Galaxies used in this plot were 
selected to have My < —18 mag and be detected in at least 4 bands. 
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Fig. 2. — a) Left : The distribution of rest frame My for galaxies in our sample is shown 
as a function of redshift. The lines indicate the expected absolute magnitudes for L* and 
L*/10 galaxies assuming passive evolution brightening of 1.2 mag from z = to 1. The 
horizontal line at My = — 18 corresponds to the absolute magnitude cutoff used here. The 
lower envelope in the grayscale is imposed by our sample cutoff at Iab < 25 mag. (For clarity 
only 20,000 randomly sampled galaxies are plotted.) b) Right : The distribution functions 
of stellar masses for galaxies in our sample with Iab < 25 mag and My < —18 are shown 
for redshift bins of width Az = 0.2. The higher noise seen in the low z mass functions is 
due to the much smaller volume and hence smaller number of galaxies sampled; the small 
falloff seen in the mass functions for z > 0.7 at M ste u ar < 5 x 1O 9 M is due to the apparent 
magnitude limit Iab < 25 mag for our sample. 
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Fig. 3. — The surface- density of galaxies (with background subtracted out) computed using 
the adaptive smoothing algorithm described in Appendix |A] is shown for slices at different 
redshifts. For z < 1.05, the width of the redshift slice is Az = 0.1; for the last three plots 
with z > 1.3, we use Az = 0.25. On each plot, the legend at the top gives the surface 
density (galaxies per deg -2 and per Mpc -2 ) corresponding to 1 contour unit. The number 
of contour units corresponding to the background surface density (S = -Bo) is also given. 
The scale bar indicates the size in co-moving Mpc. At higher z, we show only those redshifts 
with possibly significant structures. 
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Fig. 4. — The surface-density of galaxies in each structure found in the adaptive smoothing 
procedure. The galaxy densities are integrated in redshift over all connected 3-D pixels for 
each structure. A scale bar on each plot indicates 2 Mpc (comoving). The top legend gives 
the surface-density of galaxies corresponding to 1 contour unit and the mean redshift < z > 
of all galaxies within each structure is given. 
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Fig. 4. — cont'd 
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Fig. 4. — cont'd 
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Fig. 4. — cont'd 
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Fig. 5. — The galaxy over densities derived from the adaptive smoothing results, integrated 
in z from z = 0.25 to 1.05. The contour units are 5.1 x 10 3 galaxies deg -2 and the contours 
are at 1,2,3,4,5,6,7,8,10,12,14,16,18,20,22 and 24 units. 
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Fig. 6. — The relative locations of each of the LSS are shown with their centroid redshifts. 
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Fig. 7. — The radial distribution of galaxies is shown for each structure using radius deter- 
mined from the peak position of the number count distribution. The structures are ordered 
(as before) in decreasing numbers of pseudo 3-d space pixels. LSS above # 26 are not 
plotted since they have too few galaxies for a meaningful radial distribution. In all cases, 
the background galaxy distribution is removed. Many of the structures exhibit centrally 
peaked, well-behaved radial profiles, indicating candidate galaxy clusters with typical radii 
1-2 Mpc. In cases where the LSS has multiple peaks of comparable amplitude, the radial 
distribution appears non-monotonic. 
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Fig. 8. — a) Left : The distribution of total stellar masses for the structures are shown 
for masses 5 x 10 11 to 10 14 M Q . b) Right : The distribution of Richness measures for the 
the structures is shown, calculated as the surface density of galaxies brighter than 2 mag 
below the 3rd brightest galaxy (see text). Specifically, Richness is defined as the number 
of galaxies within the central R < 1 Mpc brighter than 2 magnitudes fainter than the 3rd 
brightest galaxy within the cluster (i.e. My < My (3rd) + 2 ). Radius is measured from the 
location of peak surface density as in Figure [71 
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Fig. 9. — The redshift distributions of summed structure masses are shown for redshift 
bins Az = 0.25. Uncertainties due to Poisson and cosmic variance are comparable and the 
expected total relative variance in these number distributions is 40-60% (i.e. /N ~ 0.5; 
see Section [3 .5! and Table HJ). The dotted line plots the relative comoving volumes for the 
selected redshift bins, with vertical scale arbitrarily normalized - simply to indicate that 
the overall distribution of structure mass with redshift is consistent with the expected mass 
conservation. The apparent discrepancy in the lowest redshift bin is probably not significant 
given the very large variances noted above. 
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Fig. 10. — The cumulative mass function (n(> M) = n{m) dm) is shown for the 42 
COSMOS LSS with linear mass bins of width AM* = 2 x 10 12 M . Also shown are 
the cumulative mass functio ns derived from optical and X-ray studies as summarized by 
(IReiprich fc Bohringerll2002l ). The mass function for masses within the Abell radius of each 
cluster is shown by the triangle symbols (M^); masses within the regions with density ex- 
ceeding 200p c and 500p c (where p r is the critical dens ity for the universe) are shown by the 
squares and crosses (see IReiprich fc Bohringerl 120021 ). The expected cosmic variance and 
Poisson noise are included in the error bars (see Section [3751 and Table HI). 
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Fig. 11. — The cumulative fraction of survey area with surface number density E relative to 
the average at each redshift (< E >) is shown for all redshifts and for low and high redshift 
ranges (colored curves). For comparison, the black line shows the average at z = 0.2 to 1.1 
obtained from the Millennium Simulation (see text). Jackknife tests, splitting the datasets 
in half, showed typical variances < 20%. 



-61 - 



5,0 



4,5 - 



4,0 - 



3,5 - 



3,0 - 



2,5 - 



2,0 



- 





♦ 



J J. 



I 



+ : < 1 Mpc 
a. : 1 to 5 Mpc 
q ; ove(z) 

_i I i t i 



0.2 



OA 



0,6 



1.0 



Fig. 12. — For each structure, the mean SED types of galaxies at projected R < IMpc 
and R = 1 to 5 Mpc are compared (solid line) with the mean type as a function of z (in a 
redshift bin Az = 0.2 centered on the same redshift). One set of points is shown for each 
of the structures plotted at the mean redshift of each structure. Note the very pronounced 
trend for the central Mpc in the structures to have earlier mean SED type (low type) for the 
galaxies. 
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Fig. 13. — Galaxy distributions over the areas of 6 representative LSS (ranging from the 
most massive to a quite low mass structure). The red points indicate early type (SED < 2.5) 
and the blue point later types. The larger dots are galaxies with My < —21 mag. The 
galaxy points plotted in the figures have their redshifts within the Az range given in Table 
[3] for each structure. The blank areas are either at the edge of the field or where stellar 
masking occurs. 
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Fig. 14. — The median stellar mass (M*) is shown for galaxies binned in redshift (the curves) 
as a function of environmental density. Two samples are shown : My < — 19 mag (left panel) 
and M* > 3 x 1O 9 M (right panel). (The densities have been normalized separately at each 
redshift for these plots (Equation [3]); however, above z ~ 0.25 the normalizations are similar 
and a relative density of 10 corresponds to a surface densities ~ 100 Mpc~ 2 - see text.) 
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Fig. 15. — The median absolute magnitude (My) is shown for galaxies binned in redshift 
(the curves) as a function of environmental density. Two samples are shown : My < —19 
mag (left panel) and M* > 3 x 10 9 M Q (right panel). 




Fig. 16. — The median spectral type (SED) is shown for galaxies binned in redshift (the 
curves) as a function of environmental density. Two samples are shown : My < —19 mag 
(left panel) and M,>3x 1O 9 M (right panel). 
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Fig. 17. — The median early galaxy fraction (or red galaxy fraction with type < 1.9) is 
shown for galaxies binned in redshift (the curves) as a function of environmental density. 
Two samples are shown : My < —19 mag (left panel) and M # > 3 x 1O 9 M (right panel). 
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Fig. 18. — The median star formation rate (SFR, M yr 1 ) is shown for galaxies binned 
in redshift (the curves) as a function of environmental density. Two samples are shown : 
My < -19 mag (left panel) and M t >3x 1O 9 M (right panel). 
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Fig. 19. — The median star formation timescale (t$f =M* / SFR) is shown for galaxies 
binned in redshift (the curves) as a function of environmental density. Two samples are 
shown : My < — 19 mag (left panel) and M* > 3 x 1O 9 M (right panel). 
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Fig. 20. — The median Maturity (=Ts j p/r cosmic ) is shown for galaxies binned in redshift (the 
curves) as a function of environmental density. M < 1 correspond to active, star forming 
galaxies and M > 1 corresponds to galaxies which formed most of their stars at earlier epochs 
but see text for a more detailed discussion for interpretation of the Maturity. 
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Fig. 21. — The Maturity (= TsF/Tcosmic) is shown as a function of galactic mass for high 
and low density environments in three redshift bins. The cut between high and low density 
was taken at relative density p = 45, i.e. between the middle two bins of the 4 density bins 
used earlier. M < 1 correspond to active, star forming galaxies and M > 1 corresponds to 
galaxies which formed most of their stars at earlier epochs but see text for a more detailed 
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Fig. 22. — a) Top : The projected model densities of the clusters are shown for the parameters 
given in Table The number to the upper left of each feature are the structure ID's given 
in Table [5] The contour units are given in the top legend, along with the background S 
in the same units, b) Bottom : The positions of 10,000 galaxies (including 5260 randomly 
placed galaxies) are shown for the cluster distribution in the top figure. 
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Fig. 23. — The density distribution recovered using the adaptive smoothing algorithm. 
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Fig. 24. — Application of the adaptive smoothing algorithm to the dark matter p article 
distribution in one of the Virgo consortium ACDM simulations (IBenson et al.ll200ll ). The 
algorithm detects virtually all significant structures seen visually in the simulation. More 
importantly, the characterisitcs of the detected structures (their scales and changes as a 
function of redshift) correspond closely to those detected in the galaxy overdensities from 
the photometric redshift catalog. 
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Fig. 25. — Application of the adaptive smoothing algorithm to the galaxies generated in the 
Millennium Simulation COSMOS light cone were processed identically to the observed galaxy 
distribution, including the imposition of redshift uncertainties like those in the observations 
and using the same smoothing parameters. The figure shows the projected distribution of 
overdensities integrated from z = 0.2 to 1.05 and with contours identical to those employed 
for the observations in Figure [51 



